Merging White Dwarfs and Thermonuclear 

Supernovae 

M. H. van Kerkwijk 

Department of Astronomy & Astrophysics, University of Toronto, 50 Saint George Street, 
Toronto, ON, M5S 3H4, Canada; mhvk@ astro, utoronto.ca 

Thermonuclear supernovae result when interaction with a companion reignites nuclear fusion in 
a carbon-oxygen white dwarf, causing a thermonuclear runaway, a catastrophic gain in pressure, 
and the disintegration of the whole white dwarf. It is usually thought that fusion is reignited 
in near-pycnonuclear conditions when the white dwarf approaches the Chandrasekhar mass. I 
briefly describe two long-standing problems faced by this scenario, and our suggestion that 
these supernovae instead result from mergers of carbon-oxygen white dwarfs, including those 
that produce sub-Chandrasekhar mass remnants. I then turn to possible observational tests, in 
particular those that test the absence or presence of electron captures during the burning. 
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1. The Current Paradigm and Its Problems 

A thermonuclear or type la supernova (SN la) is generally thought to be produced 
by a carbon-oxygen white dwarf that accretes matter relatively slowly, on timescales 
of > 10 6 yr (limited by the rate at which heat from accretion and pos sible nuclear 



processing can be radiated, viz., the Eddington luminosity; for reviews, iNomoto et al. 



19941 ; iHillebrandt & Niemeyerll2000h . As the white dwarf accretes, its interior is heated, 
but it does not reach ignition, because at temperatures of > 10 8 K, neutrino cooling 
becomes efficient enough to balance the heating (see Fig. [T]). However, as the white 
dwarf approaches the Chandrasekhar mass, the density in its core becomes so high 
that fusion becomes possible at lower temperatures (in partly pycno-nuclear conditions; 
Fig. [U. Once this happens, a runaway ensues, stopping only when degeneracy is lifted 
and thermal pressure can expand and cool the region. The process triggers a burning front 
that proceeds through the white dwarf, generating the energy that eventually disrupts it. 

The above is physically plausible, but it has two well-known problems. I briefly 
describe these below, before turning to our alternative. 



(a) The Paucity of Possible Progenitor Systems 

Over the age of the U niverse, for every sola r mass of stars form ed, ~ 0.0023 ±0.0006 
SN la seem to occur dMannucci et all 120051 : iMaoz et all 1201 lh . Since ~ 0.22 white 
dwarfs are expected for every solar mass formed (the remainder being in low-mass stars 
that are still alive), one infers that a surprisingly high fraction, of ~ 1%, of all white 
dwarfs eventually produce SN la. Comparing different galaxies, th e instantaneous SN l a 
rate similarly seems to be ~ 1% of the white-dwarf formation rate dPritchet et a/.L 120081) . 
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Most SN la models invoke "single dege nerate" progenitors, in which a white dwarf 
accretes from a non-degenerate companion (Whelan & IbenL Il973l) . In principle, ample 
numbers of such binaries exist and several routes to explosions have been proposed 
(Iben . Jr. & TutukovL [T984). ^o route, however, seems both common and efficient. 

The main problem is that if mass transfer is slow, unstable hydrogen fusion in 
the accreting matter c auses novae, which in mo st cases appear to remove as much 
mass as was accreted dTownsley & Bildstenl 120041 : t hough white dwarfs in cataclysmic 
variables are more massive than in their progenitors, IZorotovic et a/.ll201lh . If accretion 
is faster, hydrogen burns sta bly, but only in a sma ll range of accretion rate can expansion 
and mass loss b e avoided dNomoto et al\ [20071 : for the effect of helium flashes, see 
Idan et al\ 12012 ). Empi rically, the best-suited systems are the supersoft X- ray sources 



(Rappap ort et a/1 1 19941). but those are far too rare to explain the SN la rates (iDi Stefanol. 
2010af Gilfanov & Bogdarj, 120101) . We may be missing syste ms, e.g., more rapidl y 
accreting white dwarfs that expanded and hid from X-ray view (lHachisu et all HblO). 
However, for such sources - as for many single-degenerate channels - the lack of 
evidence for (entrained) hydrogen in SN la is surpris ing (unless the explosion can 



some how be delayed, as in the "spin-up/down" model: |Justhamll201 ll : IDi Stefano et al. 



2011). 



Another class of S N la models invoke "double degenerates," where a white dwarf 
merges with another (IWebbinkL 1 19841: llben. Jr. & TutukovL [l984). As ignition is not 
expected during the merger (except perhaps for unusually massive, > O.9M white 
dwarfs, Pak mor et al\ |2012|) . it is usually assumed an explosion will follow only if 
the combine d mass exc e eds the Chandrasekhar mass. This is rare, however, and both 
theoretical dRuiter et all 120091: iMennekens et all 1201 dl : Ivan Kerkwij k et all 20101) and 
empirical dBadenes & MaozL 1201 2l) rate estimates fall well below the SN la rate. 
Furthermore, the observed number of supersoft symbiotic proge nitors with t he requ ired 
massive white dwarfs is substantially smaller than that expected (IDi StefanoLl2010M) . 



(b) The Difficulty of Reproducing SN la Properties 

In degenerate matter, a thermonuclear runaway will proceed to completion 
unless degeneracy is lifted, and thermal pressure can expand and cool matter. 
After initial ignition, what happens depends on the conditions. For sufficiently high 
overpressure in a sufficiently larg e region (where what is "sufficient" remains to be 
understood; Seitenza hl et all 12009). a detonation is triggered: a shock strong enough 



to cause neighbouring matter to ignite and burn in turn. Since a detonation proceeds 
supersonically, the white dwarf has no time to expand and the initial density everywhere 
determines the end-point of the runaway. For a near-Chandra sekhar mass white dwarf, 
most matter is at very high density and thus far too much 56 Ni is produced. 

For a near-Chandrasekhar mass white dwarf, however, the energy release even from 
fusion up to 56 Ni does not lead to strong overpressure, and a deflagration is more likely, 
where neighbouring regions are ignited by a heat wave rather than a shock. Since a 
deflagration is sub-sonic, the white dwarf expands as the burning front progresses. Thus, 
burning takes place at lower density, reaching lower peak temperatures and producing 
less 56 Ni. Unfortunately, the burning front a ppears to be too slow, making it impossible 



to produce sufficiently energetic explosions (Hillebra ndt & Niemeyerl 12000). 



Another problem is that both pure detonation and pure deflagration models do not 
naturally reproduce the range in SN la properties, which trace a (nearly) single-parameter 
family, reflecting a roughly factor 5 range in the amount of 56 Ni that is synthesized 
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Figure 1. Temperate-density tracks leading to thermonuclear runaways. In the standard picture (on the right), 
a carbon-oxygen white dwarf accretes slowly, on a > 10 6 yr timescale, and neutrino cooling keeps its internal 
temperature below a few 10 8 K (see the contours of constant cooling time t v ). When it approaches the 
Chandrasekhar limit and its central density rises dramatically, carbon fusion is ignited (where Tec = tv. i- e -, 
the fusion heating time matches the neutrino cooling time), and - after 10 0-1000 yr of simmering - a thermal 
runaway ensues (at roughly constant pressure). In our alternative picture (van Kerkwii k et at .,2010), a white 
dwarf merger leads to a rapidly ro tating remnant with a temperature-density profile like that shown on the left 
(fromToren- Aguila r et fll.ll2009l . for a merger of two 0.6M S white dwarfs). This is initially not hot enough 
for ignition, but as the remnant disk accretes or the core spins down, the interior will be compressed and 
heated roughly adiabatically, until carbon fusion becomes faster than the accretion or spin-down timescales 
and the thermonuclear runaway starts (along a curved constant-pressure contour, as degeneracy is lifted). 



(iPhiHipsl . 119931 : iMazzali et ali [2007). The above issues can be resolved with an ad hoc 
assumption, that an initial deflagration transitions into a detonation ( Khokhlo vl. 1 1 99 lh . 
If so, the timing of the transition could determine how far the white dwarf expanded 
and thus how much 56 Ni was produced. Even with this assumption, however, it remains 
unclear why the outcome would depend on the population which the progenitor is in, i.e., 
why, as is observed, more luminous SN la preferentially occur in younger populations 
(lHamuv et a/.LIl995l : ISurlivan et a/.Ll2010h . 



2. Sub-Chandrasekhar Mass Mergers as SN la Progenitors? 

SN la could be understood more easily if they arose from sub-Chandrasekhar white 
dwarfs. Since for increasing mass, a larger fraction is dense enough to produce 56 Ni (p > 
10 7 gcm~ 3 ), a range of 56 Ni mass would be expected. Also, since more massive white 
dwarfs are the progeny of shorter-lived stars, younger populations should preferentially 
host luminous SN la. Encouragingly, pure detonations of white dwarfs with masses 
between 0.9 and 1.2M Q reproduce the range in SN la properties, including, roughly, 
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their lightcurves and spectra (Shigeva ma et all 1 19921 ; ISim et all l2010h . Not clear yet, 
however, is whether the distribution in luminosity can also be matched easily. 

The difficulty for sub-Chandrasekhar white dwarfs is to get them hot enough to 
ignite. To overcome neutrino losses, they have to be heated on a rather fast, < 10 4 yr 
timescale (see Fig. \Q. One possibility is that carbon fusion is not triggered directly, 
but indirectly, by a detonation wave started by a thermon uclear runaway in a thick 
helium layer surrounding the core (rWoosley & Weaver! Il994l) . These "double detonation" 
models, however, predict abundan ces in the outer ejecta - pro duced in the helium 
envelope - that are not seen in SN la (iHillebrandt & Niemeyerll2000l : discussions continue 
about whether these effects ca n be reduced by helium lay ers that are thinner llFink et al. 



l201d : IWooslev & Kaserj[20lTll or have mixed in carbon llKromeref aZ.ll2010ll ). Another 
possibility is that fusi on gets ignited during a merger that involves at least one massive, 



ing 

> 1M white dwarf dPakmor et all 120121) . Those, however, have expected rates even 



lower than those of super-Chandrasekhar mergers, and thus likely are too rare. 

Our alternative is that SN la result generally from merge rs of carbon-oxygen whit e 
dwarfs, including those with sub-Chan drasekhar total mass (Ivan K erkwiik ef all l2010h . 
Both theoretical (ibid.) and empirical ( Badenes & Maozl 120121) rates are a factor three 
or so higher than near-Chandrasekhar rates, making them consistent with the SN la 
rate. Furthermore, the expected range in mass matches that for which detonations yield 
sufficient Ni. The questions are whether fusion is ignited, and whether this triggers a 
detonation. 

From simulations, the outcome of a white-dwarf merger d epends strongly on whether 
the masses are similar (where "similar" is within ^0.1M ? . IZhu et aI1l201ll : Zhu et al, 
2012, in preparation). If they are not, the remnant consists of an almost unaffected core of 
the more massive white dwarf, surrounded by a hot envelope of the disrupted lower-mass 
one. For these, further evolution li kely leads to ig nition at low density, stable burning, 
and, therefore, not to a SN la (see lShen et a/.ll2012f) . 

For similar-mass white dwarfs, however, the remnants are hot throughout, and consist 
of rapidly rotating cores surrounded by thick, dense disks. Initially, the core is not hot 
enough to ignite fusion - nor dense enough to produce Ni - but as the disk accretes or 
the remnants spins down (helped by, e.g., strong magnetic fields that could be generated 
in the strongly differentially rotating remnant), it will be compressed and heated further 
(see Fig.Q]). The timescale would likely be the viscous one - hours to days - much faster 
than any relevant cooling timescale. An open question is where ignition takes place. If 
magnetic braking is important (as in a protostar or accreting pulsar), dissipation will be 
far from the remnant and ignition likely in the c ore. If accretion dominates, dissipative 
heating may lead to ignition in the outer regions (Shen et a/1.l2012l) . 



3. Observational Tests 



It seems unlikely that the question of the nature of the progenitors of SN la will be 
resolved theoretically, and hence one has to turn to observational tests. So far, most 
have focussed on trying to distinguish between the single and double-degenerate scenario, 
with conflicting results: no signature of a (former) companion in early SN la lightc urves 
(lHayden et a/.U 2010: Bia nco et all 1201 fl: iBrown et all 120121: iBloom et all 120121) or in 
SN la remnants (e.g.. ISchaefer & Pagnottall2012l: iKerzendorf et a/.ll2012l) . yet evidence 
for circumstellar medium (IPatat et all 120071 : S ternber g et all 1201 11) . 
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A different test would be to distinguish between a near or sub-Chandrasekhar mass. 
One clue is that in the near-Chandrasekhar case, where the explosion has to start with a 
deflagration, electron captures during this relatively slow phase are important, leading 
to the production o f ~0.1M Q of stable iron-peak elements, much of which is 58 Ni 
(Ma eda et a/.ll2010fll . and references therein). In contrast, for sub-Chandrasekhar models, 
where the density is much lower and the explosion has to be a fast detonation, the only 
source of the neutrons required to produce stable iron-peak elements is 22 Ne. This is 
produced during helium burning (via 14 N(a,y) 18 F(e~, v e ) I8 0(a,y) 22 Ne, where the 14 N 
is left by the CNO cycle), and wherever the temperatures become hot enough to pr oduce 
56 Ni, the excess neutrons end up mostly in 54 Fe and 58 Ni dShigeyama et qZ.L Il992h . with 
a mass of ~ (58/14)Xcno — 4% of the mass of 56 Ni; hence, the mass of 58 Ni should be 
<0.02M Q for a typical SN la with O.6M of 56 Ni. 

Given the above, an observational test would be to look for evidence for a core 
dominated by stable elements. Arguably the most direct measurement of the amount of 
58 Ni has been done from mid-infrared fine-structure lines in SN 2005df (in the nebular 
phase, when all 56 Ni has decayed; note that in other analyses often a near-Chandrasekh ar 
explosion is assumed indirectly, e.g., in using the W7 model re.g. JMazzah et a/1l2007ll ). 
These yield an estimate of ^O.OlMfn of nicke l, which is much more consistent with a 
sub-Chandrasekhar model dGerardy et Q^.ll2007l : note that these authors argued even this 
small mass was evidence for electron captures, but they did not consider the effect o f 
Ne). Similarly, the meteoritic abundance of nickel is ~ 5% that of iron dCoxL [2000), 
which is more easily understood in sub-Cha ndrasekhar models (as already noted in, e.g., 



Shig evama ef^Zll 19921 : iNomoto et alM 994). 



In c ontrast, the presence o f an inert, colder core is inferred from flat-topped line 
profiles dMotohara et all |2006). It is unclear, however, whether this col d core reflects a 
lack of heating, or rather enhanced cooling in an "infrared catastrophe" dLeloudas et all 
120091) . Evidence for an inert cor e comes also from differences in line profiles for 
lower and higher ionisation states (IMaeda et all 1201 0&I) . differences that correlate with 
other S N la properties and are plausible for delayed detonation, near-Chandrasekhar 
models (IMaeda et all 1201 Oct) . It is not yet known what to expect for sub-Chandrasekhar 
explosions, but nebular spectroscopy nevertheless seems one of the most promising ways 
of determining whether SN la result from near or sub-Chandrasekhar mass objects. 
Ideally, one would study supernovae that cover not only a range in SN la properties 
but also in host metallicity (with which 58 Ni should scale linearly for sub-Chandrasekhar 
models; for near-C handrasekhar models, the dependence is more complicated, see, e.g., 
Jackson et al.\\20l(h . 
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